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Abstract. The aim of this article is to provide a systematic analysis of 
the conditions such that Fourier transform valuation formulas are valid 
in a general framework; i.e. when the option has an arbitrary payoff 
function and depends on the path of the asset price process. An interplay 
between the conditions on the payoff function and the process arises 
naturally. We also extend these results to the multi-dimensional case, 
and discuss the calculation of Greeks by Fourier transform methods. As 
an application, we price options on the minimum of two assets in Levy 
and stochastic volatility models. 



1. Introduction 

Since the seminal work: of Carr and Madan (1999) and Raible (2000) on 
the valuation of options with Fourier transform methods, there have been 
several articles dealing with extensions and analysis of these valuation for- 
mulas. This literature focuses on the extension of the method to other situ- 
ations, e.g. the pricing of exotic or multi-asset derivatives, or on the analysis 
of the discretization error of the fast Fourier transform. 

The article of Borovkov and Novikov (2002) deals with the application of 
Fourier transform valuation formulas for the pricing of some exotic options, 
while Hubalek et al. (2006) use similar techniques for hedging purposes. Lee 
(2004) provides an analysis of the discretization error in the fast Fourier 
transform, while Lord (2008) extends the method to the pricing of options 
with early exercise features. Recently, Hubalek and Kallsen (2005), Biagini 
et al. (2008) and Hurd and Zhou (2009) extend the method to accommo- 
date options on several assets, considering basket options, spread options 
and catastrophe insurance derivatives. Dufresne et al. (2009) also consider 
the valuation of payoffs arising in insurance mathematics by Fourier meth- 
ods. In addition, the books of Cont and Tankov (2003), Boyarchenko and 
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Levendorskii (2002) and Schoutens (2003) are also discussing Fourier trans- 
form methods for option pricing. Let us point out that all these results are 
intimately related to Parseval's formula, cf. Katznelson (2004, VI. 2. 2). 

The aim of our article is to provide a systematic analysis of the condi- 
tions required for the existence of Fourier transform valuation formulas in 
a general framework: i.e. when the underlying variable can depend on the 
path of the price process and the payoff function can be discontinuous. Such 
an analysis seems to be missing in the literature. 

In their work, Carr and Madan (1999), Raible (2000) and most others are 
usually imposing a continuity assumption, either on the payoff function or 
on the random variable (i.e. existence of a Lebesgue density). However, when 
considering e.g. a one-touch option on a Levy-driven asset, both assumptions 
fail: the payoff function is clearly discontinuous, while a priori not much is 
known about the existence of a density for the distribution of the supremum 
of a Levy process. Analogous situations can also arise in higher dimensions. 

The key idea in Fourier transform methods for option pricing lies in the 
separation of the underlying process and the payoff function. We derive con- 
ditions on the moment generating function of the underlying random vari- 
able and the Fourier transform of the payoff function such that Fourier based 
valuation formulas hold true in one and several dimensions. An interesting 
interplay between the continuity conditions imposed on the payoff function 
and the random variable arises naturally. We also derive a result that allows 
to easily verify the conditions on the payoff function (cf. Lemma 2.5). 

The results of our analysis can be briefly summarized as follows: for gen- 
eral continuous payoff functions or for variables, whose distribution has a 
Lebesgue density, the valuation formulas using Fourier transforms are valid 
as Lebesgue integrals, in one and several dimensions. When the payoff func- 
tion is discontinuous and the random variable might not possess a Lebesgue 
density then, in dimension one, we get pointwise convergence of the val- 
uation formulas under additional assumptions, that are typically satisfied. 
In several dimensions pointwise convergence fails, but we can deduce the 
valuation function as an L^-limit. 

In addition, the structure of the valuation formulas allows us to derive 
easily formulas for the sensitivities of the option price with respect to the 
various parameters; otherwise, Malliavin calculus techniques or cubature for- 
mulas have to be employed, cf. e.g. Fournie et al. (1999), Teichmann (2006) 
and Kohatsu-Higa and Yasuda (2009) . We discuss results regarding the sen- 
sitivities with respect to the initial value, i.e. the delta and the gamma. It 
turns out that the trade-off between continuity conditions on the payoff 
function and the random variable established for the valuation formulas, 
becomes now a trade-off between integrability and smoothness conditions 
for the calculation of the sensitivities. 

The valuation formulas allow to compute prices of European options very 
fast, hence they allow the efficient calibration of the model to market data 
for a large variety of driving processes, such as Levy processes and affine 
stochastic volatility models. Indeed, for Levy and affine processes the mo- 
ment generating function is usually known explicitly, hence these models are 
tailor-made for Fourier transform pricing formulas. 
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We also mention here that the Fourier transform based approach can be 
apphed for the efficient computation of prices in other frameworks as well. 
An important area is the valuation of interest rate derivatives in Levy driven 
models. Levy term structure models were developed in a series of papers in 
the last ten years; this development is surveyed in Eberlein and Kluge (2007). 
For the Fourier based formulas we mention the two papers by Eberlein and 
Kluge (2006a, 2006b), where caps, floors, and swaptions as well as interest 
rate digital and range digital options are discussed; furthermore Eberlein and 
Koval (2006), where cross currency derivatives are considered and Eberlein, 
Kluge, and Schonbucher (2006), where pricing formulas for credit default 
swaptions are derived. Moreover, in the framework of the 'affine LIBOR' 
model (cf. Keller-Ressel et al. 2009) caps and swaptions can be easily priced 
by Fourier based methods. 

This paper is organized as follows: in Section 2 we present valuation for- 
mulas in the single asset case, and in Section 3 we deal with the valuation 
of options on several assets. In Section 4 we discuss sensitivities. In Section 
5 we review examples of commonly used payoff functions, in dimension one 
and in multiple dimensions. In Section 6 we review Levy and affine pro- 
cesses. Finally, in Section 7 we provide numerical examples for the valuation 
of options on several assets in Levy and affine stochastic volatility models. 

2. Option valuation: single asset 

1. Let "B = (r2,JF, F,P) be a stochastic basis in the sense of Jacod and 
Shiryaev (2003, 1.1.3), where = Tt and F = {Tt)o<t<T- We model the 
price process of a financial asset, e.g. a stock or an FX rate, as an exponential 
semimartingale S = {St)o<t<T, i-e. a stochastic process with representation 

St = Soe^\ 0<t<T (2.1) 

(shortly: S = Sqc^), where H = {Ht)o<t<T is a semimartingale with Hq = 0. 
Every semimartingale H = {Ht)o<t<T admits a canonical representation 

H = B + H^ + h{x) *{n-v) + {x- h{x)) * fi, (2.2) 

where h = h[x) is a truncation function, B = {Bt)o<t<T is a predictable 
process of bounded variation, = {H^)Q<:t<T is the continuous martingale 
part of H with predictable quadratic characteristic {H'^) = C, and v is 
the predictable compensator of the random measure of jumps fi of H. Here 
W * n denotes the integral process of W with respect to n, and W * {n — 
v) denotes the stochastic integral of W with respect to the compensated 
random measure fJ- — ly; cf. Jacod and Shiryaev (2003, Chapter II). 

Let J^{P), resp. ^Aloc{P), denote the class of all martingales, resp. local 
martingales, on the given stochastic basis "B. 

Subject to the assumption that the process l{x>i}e^ * has bounded 
variation, we can deduce the martingale condition 

S = Soe" G MiociP) 4^ 5 + ^ + (e^' - 1 - h{x)) * i/ = 0; (2.3) 

cf. Eberlein et al. (2008) for details. The martingale condition can also be 
expressed in terms of the cumulant process K associated to {B^C^v), i.e. 
K{1) = 0; for the cumulant process see Jacod and Shiryaev (2003). 
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Throughout this work, we assume that P is an (equivalent) martingale 
measure for the asset S and the martingale condition is in force; moreover, 
for simplicity we assume that the interest rate and dividend yield are zero. By 
no-arbitrage theory the price of an option on S is calculated as its discounted 
expected payoff. 

2. Let Y = (lt)o<t<r be a stochastic process on the given basis. We denote 
by y = {Yt)o<t<T and F = {Y_^)o<t<T the supremum and the infimum 
processes of Y respectively, i.e. 

Yt = sup Yu and Y_^ = inf 1^. 

0<M<t 0<«<i 

Notice that since the exponential function is monotonically increasing, 
the supremum processes of S and H are related via 

St= sup (5oe^') = 5oe""Po<'<^^' = 5oe^^. (2.4) 

0<t<T 

Similarly, the infimum processes of S and H are related via 

St = Soe^T, 

3. The aim of this work is to tackle the problem of efficient valuation for 
plain vanilla options, such as European call and put options, as well as for 
exotic path-dependent options, such as lookback and one-touch options, in a 
unified framework. Therefore, we will analyze and prove valuation formulas 
for options on an asset S = Soe^ with a payoff at maturity T that may 
depend on the whole path of S up to time T. These results, together with 
analyticity conditions on the Wiener-Hopf factors, will be used in the com- 
panion paper (Eberlein, Glau, and Papapantoleon 2009) for the pricing of 
one-touch and lookback options in Levy models. 

The following example of a fixed strike lookback option will serve as a 
guideline for our methodology; note that using (2.4) it can be re-written as 

(5t-/^)+ = (Soe^^-if)+. (2.5) 

In order to incorporate both plain vanilla options and exotic options in 
a single framework we separate the payoff function from the underlying 
process, where: 

(a) the underlying process can be the log-asset price process or the supre- 
mum/infimum of the log-asset price process or an average of the log-asset 
price process. This process will always be denoted by X (i.e. X = H or 
X = H OT X = H, etc.); 

(b) the pay ojf function is an arbitrary function / : M ^ M_|_, for example 
fix) = (e^ - K)+ or fix) = 1{,.>b}, for K,Bg R+. 

Clearly, we regard options as dependent on the underlying process X, i.e. 
on (some functional of) the logarithm of the asset price process S. The main 
advantage is that the characteristic function of X is easier to handle than 
that of (some functional of) S; for example, for a Levy process i7 = X it is 
already known in advance. 

Moreover, we consider exactly those options where we can incorporate the 
path-dependence of the option payoff into the underlying process X. Euro- 
pean vanilla options are a trivial example, as there is no path-dependence; a 
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non-trivial, example are options on the supremum, see again (2.4) and (2.5). 
Other examples are the geometric Asian option and forward-start options. 

In addition, we will assume that the initial value of the underlying process 
X is zero; this is the case in all natural examples in mathematical finance. 
The initial value 5*0 of the asset price process S plays a particular role, 
because it is convenient to consider the option price as a function of it, or 
more specifically function oi s = — log Sq. 

Hence, we express a general payoff as 

^(5oe^SO<t<r) =/(Xt-.), (2.6) 

where / is a payoff function and X is the underlying process, i.e. an adapted 
process, possibly depending on the full history of i?, with 

Xt:=^{Hs,G<s<t) fortG[0,r], 

and ^ a measurable functional. Therefore, the time-0 price of the option is 
provided by the (discounted) expected payoff, i.e. 

Nf{X-s) = E[^[SuQ<t<T)] =E[f[XT-s)]. (2.7) 

Note that we consider 'European style' options, in the sense that the 
holder or writer do not have the right to exercise or terminate the option 
before maturity. 

Remark 2.1. In case the interest rate r and the dividend yield 5 are non- 
zero, then the martingale condition (2.3) reads 

{5-r)t + Bt + ^ + (e^- - 1 - h{x)) * ut = ^ 

for all t, and the option price is given by Y f{X; s) = E\^f{XT — s)'j . 

4. The first result focuses on options with continuous payoff functions, such 
as European plain vanilla options, but also lookback options. 

Let denote the law, Mx^ the moment generating function and ^pxt 
the (extended) characteristic function of the random variable Xt', that is 

Mx^(n) = ^[e"^^] =ipxA-iu), 

for suitable u £ C For any payoff function / let g denote the dampened 
payoff function, defined via 

gix) = e-^f{x) (2.8) 

for some i? E M. Let g denote the (extended) Fourier transform of a function 

5, and LJj^(M) the space of bounded, continuous functions in L^(M). 

In order to derive a valuation formula for an option with an arbitrary 
continuous payoff function /, we will impose the following conditions. 

(CI): Assume that g G ivbc(^)- 
(C2): Assume that Mxj.{R) exists. 
(C3): Assume that g € L^(R). 

Theorem 2.2. // the asset price process is modeled as an exponential semi- 
martingale process according to (2.1)-(2.3) and conditions (C1)-(C3) are in 
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force, then the time-0 price function is given by 

Yf{X; s) = / e-™VxT {u - iR)f{iR - u)du. (2.9) 

2vr J 

R 

Proof. Using (2.7) and (2.8) we have 

Nf{X-s) = j f{XT-s)dP = e-'^' J e%(x-.)Px^(dx). (2.10) 
n R 
By assumption (CI), g G L^(]R) and the Fourier transform of g, 

g{u) = [ e»-5(x)dx, 



is well defined for every u G M and is also continuous and bounded. Ad- 
ditionally, using assumption (C3) we immediately have that g E LJj^(M). 
Therefore, using the Inversion Theorem (cf. Deitmar 2004, Theorem 3.4.4), 
g can be inverted and g can be represented, for all x G M, as 

9i^) = ^ [ e-'^-g{u)du. (2.11) 



27r 

k 

Now, returning to the valuation problem (2.10) we get that 

V/(X;.) =e-^-y" e^^l^ J e-^(^-^)"?(n)dn J Px^(dx) 

J e^"M J e*(-"-*^)^Pxr(dx) j?(tx)du 



27r 
"2^ 



j e''''ipxT{-u-iR)f{u + iR)du, (2.12) 



where for the second equality we have applied Fubini's theorem; moreover, 
for the last equality we have 

g(u) = J e™^e-^^/(x)dx = f{u + iR). 

R 

Finally, the application of Fubini's theorem is justified since 
J J e«-|e-^"(^-^)||?(n)|d^zPx^(dx) < J e^^^l^J \g{u)\du^ PxA^^) 

R R R R 

< KMxt{R) < oo, 

where we have used again that 'g G L"'^(R), and the finiteness of Mx^iR) is 
given by Assumption (C2). □ 

Remark 2.3. We could also replace assumptions (CI) and (C3) with the 
following conditions 

(CI'): g e Li(]R) and (C3'): e^^^^^ e LH 
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Condition (C3') yields that e^^Pxj, possesses a continuous bounded Lebesgue 
density, say p; cf. Breiman (1968, Theorem 8.39). Then, we can identify p, 
instead of g, with the inverse of its Fourier transform and the proof goes 
through with the obvious modifications. This statement is almost identical 
to Theorem 3.2 in Raible (2000). 

Remark 2.4 (Numerical evaluation). The option price represented as an 
integral of the form (2.9) can be evaluated numerically very fast. The fol- 
lowing simple observation can speed up the computation of this expression 
even further: notice that for a fixed maturity T, the characteristic function 
- which is the computationally expensive part - should only be evaluated 
once for all different strikes or initial values. The gain in computational 
time will be significant when considering models where the characteristic 
function is not known in closed form; e.g. in affine models where one might 
need to solve a Riccati equation to obtain the characteristic function. This 
observation has been termed 'caching' by some authors (cf. Kilin 2007) 

5. Apart from (C3), the prerequisites of Theorem 2.2 are quite easy to check 
in specific cases. In general, it is also an interesting question to know when 
the Fourier transform of an integrable function is integrable. The problem is 
well understood for smooth (C^ or C°°) functions, see e.g. Deitmar (2004), 
but the functions we are dealing with are typically not smooth. Hence, we 
will provide below an easy-to-check condition for a non-smooth function to 
have an integrable Fourier transform. 

Let us consider the Sobolev space H^{M.), with 

if^(M) = G l2(M) dg exists and dg G l2(M)|, 

where dg denotes the weak derivative of a function g; see e.g. Sauvigny 
(2006). Let g £ H^(W), then from Proposition 5.2.1 in Zimmer (1990) we 
get that 

dg{u) = —iug{u) (2-13) 

and g,dg e L'^{M.). 

Lemma 2.5. Let g G H^{R), then g G L^{R). 
Proof. Using the above results, we have that 

oo>y" (^\g{u)f + \dg{u)f^du = j \g{u)\^ [l + \u\^)du. (2.14) 

R R 

Now, by the Holder inequality, using (1 -|- |n|)^ < 3(1 -|- |up) and (2.14), we 
get that 



\g{u)\du = J |?(^)| I ^ j^j dn 



i_ 1 

'2 / r 1 \ 2 



M R 

and the result is proved. □ 
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Remark 2.6. A similar statement can be proved for functions in the Sobolev- 
Slobodeckij space ff'*(M), for s > ^. 

6. Next, we deal with the valuation formula for options whose payoff func- 
tion can be discontinuous, while at the same time the measure Pxj. does not 
necessarily possess a Lebesgue density. Such a situation arises typically when 
pricing one-touch options in purely discontinuous Levy models. Hence, we 
need to impose different conditions, and we derive the valuation formula as 
a pointwise limit by generalizing the proof of Theorem 3.2 in Raible (2000). 
A similar result (Theorem 1 in Dufresne et al. 2009) has been pointed out 
to us by one of the referees. 

In this and the following sections we will make use of the following nota- 
tion; we define the function g and the measure g as follows 

g{x) := g{-x) and g{dx) := e^Px^idx). 

Moreover g{M.) = J g{dx), while g*g denotes the convolution of the function 
g with the measure g. In this case we will use the following assumptions. 
(Dl): Assume that g e L^{M.). 

(D2): Assume that Mxj,{R) exists (^=^> g{R) < co). 

Theorem 2.7. Let the asset price process be modeled as an exponential 
semimartingale process according to (2.1) -(2.3) and conditions (D1)-(D2) 
he in force. The time-0 price function is given by 

R 

Yf{X;s) = lim / e-''''ipxT{u-iR)f{iR-u)du, (2.15) 

A— >cxD 27r J 
-A 

at the point s gM., ifYf{X;-) is of bounded variation in a neighborhood of 
s, andYf{X;-) is continuous ats. 

Remark 2.8. In Section 5 we will relate the conditions on the valuation 
function V/ to properties of the measure Pxj. for specific (dampened) payoff 
functions g. These properties are easily checkable - and typically satisfied - 
in many models. 

Proof. Starting from (2.10), we can represent the option price function as a 
convolution of g and g as follows 



Yf{X-s) = e-^^ j e"^g{x - s)PxA^x) = e"^' j g{s - x)g{dx) 

= e-'^'g*g{s). (2.16) 

Using that g G -L^(M), hence also g G L^(R), and gi^) < oo we get that 
g * g G L^(R), since 

\\9 * qWl^m.) < q(^) MIlHr) < oo; (2-17) 

compare with Young's inequality, cf. Katznelson (2004, IV. 1.6). Therefore, 
the Fourier transform of the convolution is well defined and we can deduce 
that, for all u G M, 

f*~g{u) =|(n) • g{u); 
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compare with Theorem 2.1.1 in Bochner (1955). 

By (2.17) we can apply the inversion theorem for the Fourier transform, 
cf. Satz 4.2.1 in Doetsch (1950), and get 

A 

\ {g * q{^+) + 9 * q{^-)) = 7^ lim / e-'^'Q{uMu)du, (2.18) 

-A 

if there exists a neighborhood of s where s h-> g*Q{s) is of bounded variation. 

We proceed as follows: first we show that the function ^ i-^ ^ =f: g{s) has 
bounded variation; then we show that this map is also continuous, which 
yields that the left hand side of (2.18) equals g * q{s)- 

For that purpose, we re- write (2.16) as 

then, ^ * ^> is of bounded variation on a compact interval [a, b] if and only if 
Yf{X; •) £ BV{[a, b]); this holds because the map s i— > e^^ is of bounded vari- 
ation on any bounded interval on M, and the fact that the space BV{[a, b]) 
forms an algebra; cf. Satz 91.3 in Heuser (1993). Moreover, .? is a continuity 
point of ^ * ^3 if and only if Yf{X; •) is continuous at s. 
In addition, we have that 

|(u) = J e-^'^^e-^^/(x)dx = f{iR - u) (2.19) 
R 

and 

g(u) = J e'^'^'e^Pxridx) = tpxriu - iR). (2.20) 

Hence, (2.18) together with (2.19), (2.20) and the considerations regarding 
the continuity and bounded variation properties of the value function yield 
the required result. □ 

3. Option valuation: multiple assets 

1. We would like to establish valuation formulas for options that depend on 
several assets or on multiple functionals of one asset. Typical examples of 
options on several assets are basket options and options on the minimum or 
maximum of several assets, with payoff 

where x Ay = min{x,y}. Typical examples of options on functionals of a 
single asset are barrier options, with payoff 

(^T-K) + 1|5,>B}, 
and slide-in or corridor options, with payoff' 

N 

{St - K)^ ^ 1{l<5t,<h}, 

4=1 

at maturity T, where = Tq < Ti < • • • < Tjy = T. 
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In the previous section we proved that the valuation formulas for a single 
underlying is still valid - at least as a pointwise limit, under reasonable 
additional assumptions - even if the underlying distribution does not possess 
a Lebesgue density and the payoff is discontinuous. 

In the present section we will generalize the valuation formulas to the case 
of several underlyings. Once again, if either the joint distribution possesses 
a Lebesgue density or the payoff function is continuous, the formula is valid 
as a Lebesgue integral. In case both assumptions fail, we will encounter 
situations that are apparently of harmless nature, but where the pointwise 
convergence will fail. In this case we will establish the valuation formulas 
as an L^-limit; however, with respect to numerical evaluation, a stronger 
notion of convergence would be preferable. 

Analogously to the single asset case we assume that the asset prices evolve 
as exponential semimartingales. Let the driving process be an M"'-valued 
semimartingale H = {H^, . . . , H'^)'^ and S = {S^, . . . , S"^)""" be the vector 
of asset price processes; then each component S' of S is modeled as an 
exponential semimartingale, i.e. 

Si = S'oexpHi, 0<t<T,l<i<d, (3.1) 

where H' is an M-valued semimartingale with canonical representation 

H' = H'q + B' + H'^" + h'{x) * (/i - I/) + (x' - h'{x)) * /x, (3.2) 

with h'{x) = ejh{x). The martingale condition can be given as in eq. (3.3) 
in Eberlein, Papapantoleon, and Shiryaev (2009). 

2. In the sequel, we will price options with payoff /{Xt — s) at maturity T, 
where Xt is an .T^T-measurable M'^-valued random variable, possibly depen- 
dent of the history of the d driving processes, i.e. 

XT = 'i'{Ht,0<t<T), 

where ^ is an M'^-valued measurable functional. Further / is a measurable 
function / : ^ M+, and = (.\ . . . , s'^) E R'^ with = - log 5^. 

Analogously to the single asset case, we use the dampened payoff function 

g(x) := e-<-f^'^>/(x) for x G R"^, 

and denote by g the measure defined by 

Q{dx) ■.= e^^'''^Pxridx), 

where i? G M'^ serves as a dampening coefficient. Here (•, •) denotes the 
Euclidian scalar product in R"^. The scalar product is extended to C"' as 
follows: for u,v ^ C^, set {u,v) = "^^uiVi, i.e. we do not use the Hermitian 
inner product. Moreover, Mxj. and (px^ denote the moment generating, 
resp. characteristic, function of the random vector Xt- 

To establish our results we will make use of the following assumptions. 

(Al): Assume that g £ L^(]R^). 

(A2): Assume that Mxj-iR) exists. 

(A3): Assume that q£ L^(M^). 

Remark 3.1. We can also replace Assumptions (Al) and (A3) with the 
following assumption 
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(Al'): Assume that g G JM'^) and g G L^{W^)] 

this shows again the interplay between the continuity properties of the payoff 
function and the underlying distribution. 

Theorem 3.2. // the asset price processes are modeled as exponential semi- 
martingale processes according to (3.1)-(3.2) and conditions (Al)-(A3) are 
in force, then the time-0 price function is given by 

YfiX;s) = j e-^<"'^>Mx^(i2 + m)/(ii2-n)dt/. (3.3) 

Proof. Similarly to the one-dimensional case we have that 

Yf{X;s)=e-<'^"^g*Q{.). (3.4) 

Since g G L^(M°') and q{R'^) < oo, we get that g * g & L'^{R'^); therefore 
g * q{u) = g{u)-g{u) for all u G M.'^. By assumption we know that g G L^{M.'^); 
moreover g G L°°{W^) since |^| < ||5'||Li(iRd) < oo. These considerations yield 
that f^g G L^R'^), again by using Young's inequality. 

Hence, applying the formula for the Fourier inversion, cf. Corollary 1.21 
in Stein and Weiss (1971), we conclude that 

Yf{X-s) = / e-''^'^^"^t{u)g{u)du 



g-(-R,s> 
(27r)"' 



Q-^i'^'-') Mxt{R + iu)f{iR - u)<lu, 



for a.e. s G W^. 

Moreover, if 5 i-^ Yf{X; s) is continuous, then the equality holds pointwise 
for all s G M"^. The mapping (3.4) is continuous if the mapping s ^ g * g{s) 
is continuous. Using Assumption (A3) we have that g possesses a bounded 
continuous Lebesgue density p G L^{W^); cf. Proposition 2.5 (xii) in Sato 
(1999). Then g * g = g * p and 

lim q * g(s + x) = lim / q(s-\-x — z)p(z)dz 

= / li™ + V)p{x - y)dy = g* g{s) (3.5) 

yielding the continuity of the map. Note that we have used the continuity of 
p; additionally, we can interchange integration and limit using the dominated 
convergence theorem, with majorant g{-) maxx p{x). □ 

Remark 3.3. The proof using Assumption (Al') follows analogously, with 
the obvious modifications for (3.5). 

3. Next, we consider the valuation of options on several assets when the 
payoff function is discontinuous and the driving process does not necessarily 
possess a Lebesgue density. 

The main difference to the analogous situation in dimension one is that 
the pointwise convergence of capped Fourier integrals - as is the case in Satz 
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4.2.1 in Doetsch (1950) - cannot be generalized to the multidimensional case. 
M. Pinsky gives the following astonishing example to illustrate this fact, see 
section 4.1 in Pinsky (1993); let / be the indicator function of the unit ball 
in M^, then 

/ e-^<"'->/(x)dx =l-lsmiA) + oil), (3.6) 

\x\<A 

for ^ I oo. Extrapolating the convergence results from the one-dimensional 
case to M^, we would expect pointwise convergence of the spherical sum to 
the indicator function, at least in the interior of the ball; on the contrary, 
the right hand side of (3.6) is even divergent. 

As a consequence, we only derive an L^-limit for the valuation function. 
The setting is similar to the previous sections, and we need to impose the 
following conditions. 

(Gl): Assume that g £ L^W^) n L^{R'^). 
(G2): Assume that Mx^iR) exists. 

Theorem 3.4. // the asset price process is modeled as an exponential semi- 
martingale process according to (3.1)-(3.2) and conditions (G1)-(G2) are in 
force, then the time-0 price function satisfies 

Yf{X;-) = -^^L^-lmi J e"*<"'-)v9x^(n - ii?)/(iii - n)dn. (3.7) 

[~A,A]d 

Proof. Similarly to the previous section, we have that 

Yf{X;s)=e-<'^"^g*Q{.), (3.8) 

and, for all n G M'^ 

f*~Q{u) = |(n) • g{u). (3.9) 

Now, since g G L^(]R'^) n L'^{R'^), we get that | G L'^{R'^) and ||5||L2(]Rd) = 
ll5llL2(iRd); the proofs are analogous to Theorem 9.13 in Rudin (1987). More- 
over, we have that g * g € L^(R'^), because 

\\g * ^Ili2(ffid) < giM.'^f ll5lli2(K<i) < oo. 

Therefore, since also g* g £ L^(M'^) we get that g* g £ L^(M'^) n L^iW^) 
and, analogously again to Theorem 9.13 in Rudin (1987), we get that g * g € 
L'^{W^) and \\g * gh^Rd) = \\f*~g\\ i^Rdy 

Therefore, the Fourier transform in (3.9) can be inverted and the inversion 
IS given as an L^-limit; more precisely, we have 

Wg* g-i^A\\L^(Rd) ^0 (^^oo) (3.10) 

where 

l~A,A]d 

"(2^ / ^~'^''"^f(^R-^)^Xr{u-iR)du. (3.11) 

l-A,A]d 
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Finally, (3.8) and (3.10)-(3.11) yield the option price function. □ 

Remark 3.5. The problem becomes significantly simpler when dealing with 
the product /i(-'^T)/2(^r) of a continuous payoff function fi for the vari- 
able X and a discontinuous payoff function /2 for the other variable Y, even 
in the absence of Lebesgue densities. A typical example of this situation is 
the barrier option payoff, where /i(x) = (e^ — K)~^ and f2{y) = l{ea>B}- 
Then, one can make a measure change using the (normalized) continuous 
payoff as the Radon-Nikodym derivative, apply Theorem 2.2 and then The- 
orem 2.7; this leads to pointwise convergence of the valuation function. The 
measure change argument is outlined in Borovkov and Novikov (2002) and 
Papapantoleon (2007, Theorem 3.5). 



4. Sensitivities - Greeks 

The structure of the asset price model as an exponential semimartingale, 
and the resulting structure of the option price function, allows us to easily 
derive general formulas for the sensitivities of the option price with respect 
to model parameters. In this section we will focus on the sensitivities with 
respect to the initial value, i.e. delta and gamma, while sensitivities with 
respect to other parameters can be derived analogously. 

Let us rewrite the option price function as a function of the initial value, 
using that Sq = e~' , as follows: 

Yf{X- So) = ^j S^-'^'MxAR - iu)f{^ + iR)du. (4.1) 

The delta of an option is the partial derivative of the price with respect to 
the initial value. For a generic option with payoff /, we have that 



Af{X;So) 



dYf{X;So) 

1 r d 



o / ^^-S^-'''MxAR-iu)f{u + iR)Au 

R 

= ^|sr-M.,(ii-»)|^d„. (4.2, 

M 

The gamma of an option is the partial derivative of the delta with respect 
to the initial value. For a generic option with payoff /, we get 

dAf{X;So) d^Yf{X;So) 
Tf{X- So) - 525^ 

1 f oR^2~iu MxT{R-iu)J{u + iR) 
" 2^ y ^0 (i?-l-iu)-i(i?-m)-^ ^ 

R 

In the above equations we have taken for granted that we can exchange 
integration and differentiation; however, this is the crucial step and we will 
provide sufficient conditions when we are allowed to do that. Using Satz 
IV.5.7 in Elstrodt (1999) and the elementary inequality |Im/|-F|Re/| < 2|/|, 
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we get that we can differentiate under the integral sign if there exists an 
integrable function p such that for ah u G M and all Sq > 

d 

-F{u,So) < p{u), 



where 

F{u, So) = Sl^-'^'MxriR - iu)f{u + iR). 
Now we can estimate the partial derivative of the function F: 

■^F{u, So) = |e(«-i-»)i°s^o||i2 - iuWMxAR " ^^)/(« + 
obo 

< c{l + \u\)\MxAR - iu)\\f{u + iR)\ =■■ Piu). (4.4) 
Analogously we can estimate for the second derivative of F: 



F(n,5o) <c'{l + \u\^)\MxAR-iu)\\f{u + iR)\=:p'{u). (4.5) 

Sufficient conditions for the function p in (4.4), resp. p' in (4.5), to be 
integrable are that \u\\Mxj,{R — iu)\, resp. |np|Mxy(-R — iu)\, is integrable 
and /(• + iR) is bounded; the first condition dictates in particular that the 
measure Pxj, - equivalently g - has a density of class C^, resp. C^; see 
Proposition 28.1 in Sato (1999). Alternatively, a sufficient condition is that 
the function |n||/(n+ii2)|, resp. \u\'^\ f{u+iR)\, is integrable and Mxrp{R—i-) 
is bounded, highlighting once again the interplay between the properties of 
the measure and the payoff function. 

5. Examples of payoff functions 

1. Here we list some representative examples of payoff functions used in fi- 
nance, together with their Fourier transforms and comment on whether they 
satisfy some of the required assumptions for option pricing. The calculations 
for the call option are provided explicitly and for other options they follow 
along the same lines. 

Example 5.1 (Call and put option). The payoff of the standard call option 
with strike K G M+ is f{x) = (e^ - K) + . Let z G C with G (1, oo), then 
the Fourier transform of the payoff function of the call option is 

oo oo 

f{z) = y"e^"^(e^-i^)+dx= J e^^+'^^'^dx - K J e^^Mx 

R InK InK 

1 K fC^+iz 

= +K''— = — -. (5.1) 

1 + iz iz iz[l + iz) 

Now, regarding the dampened payoff function of the call option, we easily 
get for ii G (1, oo) that g G -Z^bc(^) ^ L'^W- The weak derivative of g is 

Jo, iix <lnK, 

Again, we have that dg G L^{R). Therefore, g G H^{R) and using Lemma 
2.5 we can conclude that g G L^(R). Summarizing, conditions (CI) and (C3) 
of Theorem 2.2 are fulfilled for the payoff function of the call option. 
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Similarly, for a put option, where f{x) = (K — e^)"^, we have that 

Analogously to the case of the call option, we can conclude for the dampened 
payoff function of the put option that g G LJj^(R) and g G H^(R) for R < 0, 
yielding g G L^(M). Hence, conditions (CI) and (C3) of Theorem 2.2 are 
also fulfilled for the payoff function of the put option. 

Example 5.2 (Digital option). The payoff of a digital call option with 
barrier B £ R_|_ is l{exyB}- Let z £ C with "^z G (0,oo), then the Fourier 
transform of the payoff function of the digital call option is 

/(.)=--. (5.4) 
Similarly, for a digital put option, where f{x) = Ijei^^^}, we have that 

f{z) = ^, 9zG(-cx),0). (5.5) 

For the dampened payoff function of the digital call and put option, we can 
easily check that g G L^{R) for R G (0,oo) and R G (-oo,0). 

Regarding the continuity and bounded variation properties of the value 
function, we have that 

YfiX,s)=E[l^^xr-.^j,y] =P{XT>log{B)+s) = l-Fx^{log{B)+s), 

where Fxj- denotes the cumulative distribution function of Xt- Therefore, 
s Yf{X,s) is monotonically decreasing, hence it has locally bounded 
variation. Moreover, we can conclude that s i-^ Yf{X, s) is continuous if the 
measure Px^ is atomless. 

Summarizing, condition (Dl) is always satisfied for the payoff function of 
the digital option, while the prerequisites of Theorem 2.7 on continuity and 
bounded variation are satisfied if the measure Px^ does not have atoms. 

Example 5.3. A variant of the digital option is the so-called asset-or- 
nothing digital, where the option holder receives one unit of the asset, instead 
of currency, depending on whether the underlying reaches some barrier or 
not. The payoff of the asset-or-nothing digital call option with barrier B G 
M4. is f{x) = e^l{e^>B}5 and the Fourier transform, for z € C with Qz G 
(1,00), is 

m = (5.6) 

Arguing analogously to the previous example, we can deduce that condition 
(Dl) is always satisfied for the payoff function of the asset-or-nothing digital 
option, while the prerequisites of Theorem 2.7 are satisfied if the measure 
does not have atoms. 

Example 5.4 (Double digital option). The payoff of the double digital call 
option with barriers B_,B > Q\s l{B<e^<"B}- Let z G C\{0}, then the Fourier 
transform of the payoff function is 

/(z) = ^(s^^-^-). (5.7) 
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The dampened payoff function of the double digital option satisfies g S 
L\R) for all R £ R. 

Moreover, we can decompose the value function of the double digital op- 
tion as 

Yf{X,s)=Yf,{X,s)-Yf,{X,s), 

where fi{x) = ^^^x^-^y and f2{x) = l|^<e^}. Hence, by the results of Exam- 
ple 5.2, we get that condition (Dl) is always satisfied for the payoff function 
of the double digital option, while the prerequisites of Theorem 2.7 are sat- 
isfied if the measure Px^ does not have atoms. 

Example 5.5 (Self-quanto and power options). The payoff of a self-quanto 
call option with strike K G M_|_ is f{x) = {e^ — K)~^ . The Fourier transform 
of the payoff function of the self-quanto call option, for z £ C with Q'z E 
(2,oo), is 

f{z) = -. (5.8) 

■'^ ' (l + iz)(2 + zz) ^ ' 

The payoff of a power call option with strike K G M+ and power 2 is f{x) = 
[(e^ — i^)"*"]^; for z G C with '^z G (2, oo), the Fourier transform is 

iz\ \ -I- iz)\2 -I- iz) 

The payoff functions for the respective put options are defined in the obvi- 
ous way, while the Fourier transforms are identical, with the range for the 
imaginary part of z being respectively (— cxd, 1) and (— cxo,0). 

Analogously to Example 5.1, we can deduce that conditions (CI) and 
(C3) of Theorem 2.2 are fulfilled for the payoff function of the self-quanto 
and the power option. 

Remark 5.6. For power options of higher order we refer to Raible (2000, 
Chapter 3). 

2. Next we present some examples of payoff functions for options on several 
assets and for options on multiple functionals of one asset, together with 
their corresponding Fourier transforms. 

Example 5.7 (Option on the minimum/maximum). The payoff function of 
a call option on the minimum of d assets is 

/(x) = (e^M--- Ae^'' - i^)+, 
for X G M*^. The Fourier transform of this payoff function is 

f{A = 7, 1 , (5.10) 

where z £ with Qz^ > for 1 < k < d and Q{Ylt=i ^k) > li for more 
details we refer to Appendix A. Then, we can easily deduce for the dampened 
payoff function that g G L^^^CR'^). 

Moreover, for the put option on the maximum of d assets, the payoff 
function is 



/(x) = (K - e^^i V • • • V e^ 
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for X S R'^, where aV b = max{a, b}. The Fourier transform is 

with the restriction now being Qz^ < for all 1 < A; < d. Again, we can easily 
deduce that the dampened payoff function satisfies g E Ljj^(M'^). Therefore, 
condition (Al) of Theorem 3.2 is satisfied. 

Example 5.8. A natural example of multi- asset payoff functions are prod- 
ucts of single asset payoff functions. These payoff functions have the form 

d 

/(^)=n/'(^-)' 

for X G M'^, where X; G M and /; : M ^ M+, for all 1 < t < for example, 
one can consider fi{xi) = (e^^ — K)'^ and f2{x2) = l{e^2>s}- 

The Fourier transform of these payoff functions is simply the product of 
the Fourier transform of the 'marginal' payoff functions, since 

„ d d „ d 

f(z) = J e''^^'-^l[Mx,)dx = 11 J e----/,(^.)dx, =llMz,), 

for z G C"^ and G C, with in an appropriate range such that 'g ^ [W^). 
This range, as well as other properties of /, are dictated by the corresponding 
properties of the Fourier transforms /; of the marginal payoff functions 

Remark 5.9. Further examples of multiple asset payoff functions, such as 
basket and spread options, and their Fourier transforms can be found in 
Hubalek and Kallsen (2005). 

3. We add a short remark on the rate of decay of the Fourier transform of the 
various payoff functions and its consequence for numerical implementations. 

Consider the standard call option, where the Fourier transform of the 
dampened payoff function has the form, cf. (5.1), 

j^l—R^iulogK 

9{u) = — . — , u G M. 

(it — m)[R — I — lu) 

Then, we have that 



which shows that g{u) behaves like for \u\ > 1. On the other hand, a 
similar calculation for the digital option shows that the Fourier transform 
of the dampened digital payoff behaves like - for \u\ > 1. 

Therefore, splitting a call option into the difference of an asset-or-nothing 
digital and a digital option, as many authors have proposed in the literature 
(cf. e.g. Heston 1993), is not only 'conceptually' sub-optimal, as can be seen 
by Theorems 2.2 and 2.7. More importantly, it is also not optimal from the 
numerical perspective, since the rate of decay for the digital option is much 
slower than for the call option, leading to slower numerical evaluation of the 
corresponding option prices. 



18 



E. EBERLEIN, K. GLAU, AND A. PAPAPANTOLEON 



Indeed, we have calculated the prices of call options corresponding to 11 
strikes and 10 maturities, first using the formula for the call option, and 
then representing the call option as the difference of two digital options. 
The numerical calculation using the second method lasts twice as long (6 
sees compared to less than 3 sees) in a standard Matlab implementation. 

6. Examples of driving processes 

The application of Fourier transform valuation formulas in practice re- 
quires the explicit knowledge of the moment generating function of the 
underlying random variable. As such, Fourier methods are tailor-made for 
pricing European options in Levy and affine models, since in these models 
one typically knows the moment generating function explicitly (at least up 
to the solution of a Riccati equation). In order to give a flavor, we present 
here an overview of Levy and affine processes, referring to the literature for 
specific formulas and proofs. 

In Levy processes, the moment generating function of the random variable 
is described by the celebrated Levy-Khintchine formula; for a Levy process 
H = {Ht)o<t<T with triplet (6, c. A) we have: 

£;[e<«.-f^*)] =exp(K(n) -t), (6.1) 

for suitable u £ W^, where the cumulant generating function is 

Kiu) = {b, u) + ^{u, cu) + J (e<"'"> -l-{u, h{x))^ A(dx); (6.2) 

here h denotes a suitable truncation function. The most popular Levy mod- 
els are the VG and CGMY processes (cf. Madan and Seneta 1990, Carr et 
al. 2002), the hyperbolic, NIG and GH processes (cf. Eberlein and Keller 
1995, Barndorff-Nielsen 1998, Eberlein 2001), and the Meixner model (cf. 
Schoutens and Teugels 1998). 

In affine processes, the moment generating functions are described by the 
very definition of these processes. Let X = {Xt)o<t<T be an affine process on 
the state space D = M"* x M" C M*^, starting from x ^ D; i.e., under suitable 
conditions, there exist functions (/> : [0, T] x X ^ R and tp : [0, T] x X ^ M'^ 
such that 

[e<"'^'>] = exp (cPtiu) + (M^), x)) , (6.3) 

for all (t, n, x) G [0, T] x X x D, X C M*^. The functions (p and -0 satisfy 
generalized Riccati equations, while their time derivatives 

^(■") = ^lt=o'^*('") ^(u) = -Ql\t=o'^t{u), 

are of Levy-Khintchine form (6.2); we refer to Duffie et al. (2003) and Keller- 
Ressel (2008) for comprehensive expositions and the necessary details. The 
class of affine processes contains as special cases - among others - many sto- 
chastic volatility models, such as the Heston (1993) model, the BNS model 
(cf. Barndorff-Nielsen and Shephard 2001, Nicolato and Venardos 2003), and 
time-changed Levy models (cf. Garr et al. 2003, Kallsen 2006). 
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7. Numerical illustration 

As an illustration of the applicability of Fourier-based valuation formulas 
even for the valuation of options on several assets, we present a numerical 
example on the pricing of an option on the minimum of two assets. As driving 
motions we consider a 2d normal inverse Gaussian (NIG) Levy process and 
a 2d affine stochastic volatility model. 

Let H denote a 2d NIG random variable, i.e. 

H = {H\H^) ~ NIG2(a, /3, 5, fi, A), 

where the parameters have the following domain of definition: a,5 £ ]R_|_, 
P, fi E M^, and A E M^^^ is a symmetric, positive-definite matrix; w.l.o.g. we 
can assume that det(A) = 1; in addition, > A/3). Then, the moment 
generating function of H, for it E with — (/3 + tt, A(/3 + ?x)) > 0, is 

Mh{u) = exp (^{u, fi) + 5 ( - (/?, A/3) - ^^^2 - (/3 + n, A(/? + n)))) . 

(7.1) 

In the NIG2 model, we specify the parameters a, (3,6 and A, and the drift 
vector /X is determined by the martingale condition. Note that the marginals 
H' are also NIG distributed (cf. Blassild 1981, Theorem 1), hence the drift 
vector can be easily evaluated from the cumulant of the univariate NIG 
law. The covariance matrix corresponding to the NIG2-distributed random 
variable H is 

Snig = d{a^- (/3, A/3)) (a + (a^ _ A/3)) A/3/3Ta) , 

cf. Prause (1999, eq. (4.15)). A comprehensive exposition of the multivariate 
generalized hyperbolic distributions can be found in Blassild (1981); cf. also 
Prause (1999). 

We will also consider the following affine stochastic volatility model intro- 
duced by Dempster and Hong (2002), that extends the framework of Heston 
(1993) to the multi-asset case. Let H = {H^,H^) denote the logarithm of 
the asset price processes S = (5^,5^), i.e. = log 5'; then, H', i = 1,2 
satisfy the following SDEs: 

dHj = -^afvtdt + ai^tdWl 

dHf = -^a^vtdt + a2^t^Wl 
dvt = - vt)dt + cJsi/u^dWf , 

with initial values Hq,Hq, vq > 0. The parameters have the following domain 
of definition: ai,a2,cr3 > and /i, k > 0. Here W = {W^,W'^,W^) denotes 
a 3-dimensional Brownian motion with correlation coefficients 



{W\ W^) = pi2, {W\ W'') = P13, and {W^ W^) = ^23- 
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The moment generating function of the vector H = {H^,H^) has been 
calculated by Dempster and Hong (2002); for u = {ui,U2) G we have 

/ 2C(1 - e~^*) 

MhM = exp l^{n, Ho) + - (9 - - e-<^^) ' "° 

where C = C{u), 7 = j{u), and 6 = 9{u) are 

C = ^ {oiu\ + alul + 2/>i2CJlCr2Uiti2 - CTiUl - a2'U2) , 
7 = K - - /9230-20-3'"2, 



2 -log ( ^ ^)+(^-7)t 



72 - 2aiC. 



We can deduce that all three models satisfy conditions (A2) and (A3) of 
Theorem 3.2 for certain values of R. Explicit calculations for the 2d NIG 
model are deferred to Appendix B; analogous calculations yield the results 
for the other models. 

The Fourier transform of the payoff function f{x) = (e^^ A e^^ _ K)^ , 
X S M^, corresponding to the option on the minimum of two assets is given 
by (5.10) for d = 2, and we get that condition (Al) of Theorem 3.2 is 
satisfied for i?2 > such that Ri + R2 > 1. 

Therefore, applying Theorem 3.2, the price of an option on the minimum 
of two assets is given by 

MTATiS\S^;K) = ^ J {SI)^^+''''{SI)^^+'''^MhARi + iui,R2 + iu2) 

^ dii 

(i?i + iui){R2 + iu2){Ri + R2-l + iui + iu2) 

where Mh^ denotes the moment generating function of the random vector 
Ht, and Ri,R2 are suitably chosen. 

In the numerical illustrations, we consider the following parameters: strikes 

K = {85, 90, 92.5, 95, 97.5, 100, 102.5, 105, 107.5, 110, 115} 

and times to maturity 

T= {-^,^,0.25,0.50,0.75,1.00} . 

In the 2d NIG model, we consider some typical parameters, e.g. Sq = 100, 
52 ^ 95 Q, ^ g 20, = -3.80, /32 = -2.50 and <5 = 0.150; we consider two 
matrices A"*" = ( q 1 ) and A~ = [\~2)^ which give positive and negative 
correlations respectively; indeed we get that 

^, ,0.0646 0.0191\ , ^_ / 0.0287 -0.0258\ 



NIG \ 0.0191 0.0481 J NIG 1^-0.0258 0.0556 J ' 

The option prices in these two cases are exhibited in Figure 1. 

Finally, in the stochastic volatility model we consider the parameters used 
in Dempster and Hong (2002), that is = 96, = 100, cJi = 0.5, ^2 = 1.0, 
(73 = 0.05, pi2 = 0.5, pi3 = 0.25, P23 = -0.5, vq = 0.04, k = 1.0 and 
/i = 0.04; the option prices are shown in Figure 2. 
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Figure 1 . Option prices in the 2d NIG model with positive 
(left) and negative (right) correlation. 




Figure 2. Option prices in the 2d stochastic volatility model. 

Appendix A. Fourier transforms of multi-asset options 

In this appendix we outline the derivation of the Fourier transform cor- 
responding to the payoff function of an option on the minimum of several 
assets; the derivation for the maximum is completely analogous and there- 
fore omitted. 

The payoff of a (call) option on the minimum of d assets is 

{S^ ^S'^ - K)+. 

The payoff function / corresponding to this option is given, for x € M*^, by 

f{x) = (e^'i A e^'2 A • • • A - K)+ = (e^'i^^^A-Ax, _ ^)+^ 
The following decomposition holds, if Xi ^ xj for i^j,l<i,j<d 

d d d 

i=l i=l j=l 
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where k = \ogK. Define also the auxihary functions /i, 1 < i < d, where 



fi{x) = (e^' -i^) J]l{fc<^,<^^.}. 



The dampened payoff function is g{x) = e where R G W^; 

we define analogously the dampened /j-functions, i.e. gi{x) = e~^^'^^ fi{x). 
For simplicity, we first calculate the Fourier transform of the dampened 
/i-function; for u G M'^ we get 



^ d 

Uu) = j e<»-«'->(e-^ -i^)ni{;t<x,<.,}dx 

oo oo oo 

= j j ... j e<*"-^'^> (e^^ - K)dxd . . . dxi 

k xi xi 

= / e^i^'-^^^'^^ie'^^ -K) Yl / e(^"^-^^>'^dxj- dxi 



oo 



i=2 

oo 



1 



11^=2 (-Rj { 



1 / 7^i+E-=i(i«i--Ri) j^i+E,=i(i«i--R.)' 

+ 



ni=2 i^j -iuj) \ 1 + E j=i (^""i - Rj ) Ej=i (^^j -Rj) J 



subject to the conditions Rj > for all j > 2 and Ej=i Rj > 1- 
Hence, in general we have that 



nUiRj - i^^) X 1 + E"=i(i"j- - X E"=i(^". - Rj. 



subject to the conditions Rj > for all 1 < j < (i and E^=i Rj > 1- 
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Now, we recall that /(x) = '^1=1 fiix), hence g{x) = Yli=i9ii^) which 



yields g{u) = Y.'i=i9i{u); therefore 



9[u) 



_^i+E-=i(».-fij) 



This we can also rewrite as 

f{z) = y ^, (A.l) 

subject to the conditions Qzj > for all 1 < j < d and Yl'j=i '^^j ^ ^• 

Appendix B. Calculations for the 2d NIG model 

By the moment generating function of the 2d NIG process, cf. (7.1), it is 
evident that assumption (A2) is satisfied for R £M? with — + A(/3 + 
R)) > 0. In order to verify condition (A3) we have to show that the function 
u I— > Mh{R + iu) is integrable; it suffices to show that the real part of the 
exponent of Mh{R + iu) decays like — We have 



log {Mh{R + iu)) = u) + {fi, R) + 5Va^ - {(3,A(3) 



- 5y/a'^ - {(3 + R + iu, A(/? + R + iu)) . 

Recall that the product (•, •) over is defined as follows: for u,v set 
{u,v) = Yyi'^iVi- Then 

{/3 + R + iu,A{(3 + R + iu)) 

= {j3 + R, A{P + R)) - {u, Au) + 2i{(5 + R, Au) 



and since ^/z = \J\{\z\ + '^{z)) + ij^^y^'^d-^l ~ 3?(z)), we get 
5R(log {MH{R + iu))) 

= {^i,R) +5^/a-^ - (/3,A/3) " ^{l"^ " {P + R + iu, A{P + R + iu))\ 



+ a^-{l5 + R, A{(5 + R)) + {u, Au)} 



1/2 



< {^L, R) + 5^/a^^{^^Ap) - b^o? - {[3 + R, A{[5 + R)) + (n. An) 



< {^L, R) + 5v'a'-(/5,A/3) - 6^/KZ\u\ , 
where Amin denotes the smallest eigenvalue of the matrix A. 
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